Momentum dependence of drag coefficients and heavy flavour suppression in quark 

gluon plasma 



Surasree Mazumder* , Trambak Bhattacharyya, Jan-e Alam and Santosh K Das 
Theory Division, Variable Energy Cyclotron Centre, 1/AF, Bidhan Nagar , Kolkata - 700064 

(Dated: October 10, 2011) 

The momentum dependence of the drag coefficient of heavy quarks propagating through quark 
gluon plasma (QGP) has been evaluated. The results have been used to estimate the nuclear 
suppression factor of charm and bottom quarks in QGP. We observe that the momentum dependence 
of the transport coefficients plays crucial role in the suppression of the heavy quarks and consequently 
in discerning the properties of QGP using heavy flavours as a probe. We show that the large 
suppression of the heavy quarks observed at RHIC and LHC is predominantly due to the radiative 
losses. The suppression of D° in Pb+Pb collisions at LHC energy - recently measured by the ALICE 
collaboration has also been studied. 
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I. INTRODUCTION 

Simulations of QCD equation of state (EoS) on lat- 
tice show that at very high temperatures and/or densi- 
ties the nuclear matter undergoes a phase transition to a 
new state of matter called Quark-Gluon Plasma(QGP). 
It is expected that QGP can be produced experimentally 
by colliding two nuclei at ultra-relativistic energies. Rel- 
ativistic Heavy Ion Collidcr(RHIC) at BNL and Large 
Hadron Collider(LHC) at CERN are two such experi- 
mental facilities. The depletion of hadrons with high 
transverse momentum (pr) produced in Nucleus + Nu- 
cleus collisions with respect to those produced in proton 
+ proton (pp) collisions has been considered as a signa- 
ture of QGP formation. The two main processes which 
cause the depletion are (i) the elastic collisions and (ii) 
the radiative loss or the inelastic collisions of the high en- 
ergy partons with the quarks, anti-quarks and gluons in 
the thermal bath. The importance of elastic energy loss 
in QGP diagnosis was pointed out first by Bjorken [lj. 
Later the calculations of elastic loss had been performed 
with improved techniques [3J3I and its importance was 
highlighted subsequently d, Q ■ 

The abundance of charm and bottom quarks in the 
partonic plasma, in the expected range of temperature to 
be attained in the experiments, is small. Consequently, 
the bulk properties of the plasma is not decided by them 
and hence heavy quarks may act as an efficient probe 
for the diagnosis of QGP. The collisional energy loss of 
heavy quarks Q has gained importance recently in view 
of the measured nuclear suppression in the px spectra of 
non-photonic single electrons 0, H| ■ In the present work 
we focus on the energy loss of heavy quarks in QGP in 
deducing the properties of the medium. 

Several ingredients like inclusions of non-perturbative 
contributions from the quasi- hadronic bound state [9j, 
3-body scattering effects [lfj, the dissociation of heavy 
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mesons due to its interaction with the partons in the 
thermal medium [llj and employment of running cou- 
pling constants and realistic Debye mass [HI have been 
proposed to improve the description of the experimen- 
tal data. For mass dependence of energy loss due to 
radiative processes Dokshitzer and Kharzeev |13j argue 
that radiative energy loss of heavy quarks will be sup- 
pressed compared to that of light quarks due to dead 
cone effects [14]]. However, Aurenche and Zakharov claim 
that the radiative process has an anomalous mass depen- 
dence [h| due to the finite size of the QGP which leads 
to small difference in energy loss between a heavy and 
a light quarks. Although the authors of [lfj concluded 
that the suppression of radiative loss for heavy quarks is 
due to dead cone effects but it will be fair to state that 
the issue is not settled yet. 

The other mechanism that can affect the radiative loss 
is the LPM effect [13] which depends on the relative mag- 
nitude of two time scales of the system [l8|] : the formation 
time (tf) and the mean scattering time scale (r sc ) of the 
emitted gluons. If tf > t sc then LPM suppression will 
be effective. The LPM effect is built-in in the expression 
for radiative energy loss of heavy quarks derived in [l9T - 
[22| . In the present work the LPM and the dead cone 
effects are explicitly taken into account for heavy quark 
dissipation. 

We assume here that the light quarks and gluons ther- 
malize before heavy quarks. The charm and bottom 
quarks execute Brownian motion [23l - [33j in the heat 
bath of QGP. Therefore, the interaction of the heavy 
quarks with QGP may be treated as the interactions be- 
tween equilibrium and non-equilibrium degrees of free- 
dom. The Fokker-Planck (FP) equation provide an ap- 
propriate framework for the evolution of the heavy quark 
in the expanding QGP heat bath. 

The system under study has two components. The 
equilibrium components, the QGP, is assumed to be 
formed at a temperature at an initial time Tj after 
the nuclear collisions. The QGP, due to its high in- 
ternal pressure expands, as a consequence it cools and 
reverts to hadronic phase at a temperature, T c . The 
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non-equilibrium component, the heavy quarks produced 
due to the collision of partons of the colliding nuclei has 
momentum distribution determined by the perturbative 
QCD (pQCD), which evolves due to their interaction 
with the expanding QGP background. The evolution of 
the heavy quark momentum distribution is governed by 
the FP equation. The interaction of the heavy quarks 
with the QGP is contained in the drag and diffusion co- 
efficients. The drag and diffusion coefficients are pro- 
vided as inputs, which are, in general, dependent on both 
temperature and momentum. The evolution of the tem- 
perature of the background QGP with time is governed 
by relativistic hydrodynamics. The solution of the FP 
equation at the (phase) transition point for the charm 
and bottom quarks gives the (quenched) momentum dis- 
tribution of hadrons (B and D mesons) through fragmen- 
tation. The fragmentation of the initial momentum dis- 
tribution of the heavy quarks results in the unquenched 
momentum distribution of the B and D mesons. The 
ratio of the quenched to the unquenched pt distribution 
is the nuclear suppression factor which is experimentally 
measured. The quenching is due to the dragging of the 
heavy quark by QGP. Hence the properties of the QGP 
can be extracted from the suppression factor. In contrast 
to earlier works where the momentum dependence of the 
drag coefficient were ignored (or considered its value at 
low momentum) in the present work we emphasize on its 
momentum dependence. 

The paper is organized as follows: In the next sec- 
tion we briefly describe the FP equation, in section III 
the drag coefficients for collisional and radiative processes 
have been discussed. In section IV the initial conditions 
and the space time evolution of the system have been 
described. Results have been presented in section V and 
section VI has been dedicated to the summary and dis- 
cussions. 



II. THE FOKKER PLANCK EQUATION 

The Boltzmann transport equation describing non- 
equilibrium statistical system reads: 
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where F represents external forces acting on heavy quark, 
p and E denote three momentum and energy of the probe 
respectively. In the absence of any external force in a 
uniform plasma and defining 



/( P ,t) = i j d 3 xf(x,p,t) 
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we have for the normalized probability distribution in 
momentum space, 
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For the problem under consideration we need to evalu- 
ate the collision integral for a situation where one of the 
colliding partner is in thermal equilibrium. 

To proceed in this direction, we define the collision rate 
of the heavy quark with the thermal gluons as: 



(p, k) = 7 3 



d 3 q 
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where f g is the thermal gluon distribution, fq,p is the 
relative velocity between the two collision partners, "f g is 
the degeneracy of gluon, and a 9 is the cross section for 
the heavy quark-gluon elastic interaction. The collision 
rate of the heavy quarks with the light quarks and anti- 
quarks can be evaluated in a way similar to Eq. |4] In 
terms of the transition rates the collision integral of the 
Boltzmann transport equation can written as: 

3 kHp+k, k)/(p+k)-«7(p, k)/(p)]. 

Using Landau approximation i.e. by expanding w(p + 
k, k) in powers of k and keeping upto quadratic term, 
the Boltzmann transport equation can be written as 
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where the kernels are defined as 



Ai= d 3 kw(p,k)ki 



and 



d 3 kiu(p, k)kikj 



(6) 



(7) 



(8) 



Eq. [6] is a nonlinear integro-differential equation known 
as the Landau kinetic equation. For the problem under 
consideration one of the colliding partner is in equilib- 
rium. In such situation the distribution function which 
appears in w can be replaced by thermal distribution. As 
a consequence Eq(6] becomes a linear partial differential 
equation, known as Fokker-Planck (FP) equation. As- 
suming Ai = Pijip) and Bij = D(p)5ij, the FP equation 
can be written as: 
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where the momentum, p = (p T ,p z ) = (p x ,P y ,Pz)- We 
numerically solve Eq. [9] |H[ with the boundary condi- 
tions: f(p x ,Py>t) for p x ,Py oo and the initial 
(at time t = t%) momentum distribution of charm and 
bottom quarks are taken MNR code 35]. It is evident 
from Eq. [5] that with the momentum dependence trans- 
port coefficients the FP equation becomes complicated. 
It is possible to write down the solution of the FP equa- 
tion in closed analytical form [27[ in the special case of 
momentum independent drag and diffusion coefficients. 
The computer code used for the solution of Eq. [§] has 
been used to reproduce the closed form analytical solu- 
tion of Ref. US- 



III. DRAG COEFFICIENT 

A. Collisional process 

The drag coefficient, jcoii due to collisional process can 
be written as 1261: 
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where p' = p — k and q' = q + k. The scattering matrix 
elements are given explicitly in Ref. 36]. The integra- 
tions in Eq [TBI ha s been performed using the standard 
techniques [26l l32j]. 



B. Radiative process 

The drag coefficient due to the radiative process, "f ra d 
can be related to the energy loss as follows: 



dE 
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= Irad P , 
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where p is the momentum of the particle. We evaluate 
the radiative energy loss by using the following gluon 
spectrum 371 ]: 
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where k = (fco,fcj_,fc 3 ) is the four momentum, r] = 
l/21n[(fco + fa)/ {fa — fa)] is the rapidity of the emitted 



gluon and q = (qo,q±,q3) is the four momentum of the 
exchanged (propagator) gluon. C 'a = 3 is the Casimir in- 
variant of the SU{3) adjoint representation, a s = g 2 /47r 
is the strong coupling and mp is the Debye mass. 

The dead-cone effect is taken into account through the 
factor, F [HG3 
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where #o = MJE. The average energy loss per collision, 
e is given by [H, HH 



(n„ko) 



drj d kj 
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x k 6(t sc -Tf) <d{E- k ± coshri)F 2 , (20) 

where the formation time of the emitted gluon [4p| j . 
tf = (C A /2C 2 ) 2coshry/fc_L, and C A /2C 2 = N 2 /(N 2 -1) 
for quarks with C 2 — Cf — 4/3. We perform the inte- 
grations in Eq. [20] and substitute the value of qj_ by its 
average value which is obtained as, 
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where, 
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and s fa 6ET, is the centre of mass energy squared for 
the scattering of a heavy quark with energy E off the 
thermal partons at temperature T. 

The LPM effect is taken into account by introducing 
non-zero formation time of the emitted gluon through 
the first theta function in Eq. [201 The scattering time 
scale, t sc can be estimated from the scattering rate (A) 
by using the relation, r sc = A -1 , where A is given by: 
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Note that deletion of k{ from the expression for Ai (Eq.[7]) 
gives the scattering rate, A. Knowing e and A the energy 
loss of heavy quark can be expressed as: 



dE 
di- 



ke 



(24) 



The drag due to radiative loss can now be estimated us- 
ing Eqs. [20l and l24l The effective drag due to collisional 
and radiative processes is obtained as: 7 = 7 co h + j r ad- 
It should be mentioned that the radiative process is not 
fully independent of the collisional one. Therefore, the 
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collisional and radiative transport coefficient may not be 
added to obtain the effective one. In absence of any rig- 
orous method, however, we add them up. Since the ra- 
diative loss is larger than the collisional one, therefore, 
this approximation may not be treated as unreasonable. 

A comment on the effects of Bose enhancement (BE) 
of the gluons and Pauli blocking (PB) of the quarks & 
antiquarks in the final channels on the drag is in order 
here. We have found that the average change due to 
BE and PB effects is about ~ 6% for p < 5 GeV, at 
higher p the change is negligible for the temperature do- 
main under consideration. The effects of BE and PB is 
inconsequential on the nuclear suppression factor eval- 
uated later. The diffusion coefficient has been obtained 
from the drag coefficients by using the Einstein's relation, 
D — jMT, where M is the mass of the heavy quark. It is 
important to mention at this point that the Einstein's re- 
lation is valid in the non-relativistic domain, use of this 
relation leads to underestimation of the value of D in 
the relativistic realm (23[. In our calculation the tem- 
perature dependence of the a s (4l| has been taken into 
account. 
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FIG. 1: Drag coefficients of charm assuming running strong 
coupling, as(T) and temperature dependent Debye screen- 
ing mass, m_o(r) due to its interaction with thermal gluons, 
quarks, and antiquarks. Inset: The variation of drag with p 
in the low domain. 



IV. INITIAL CONDITIONS 



In order to solve Eq. [5] we supply the initial dis- 
tribution functions, fin{PT,t) for charm and bottom 
quarks from the well known MNR code (35[. The initial 
temperature, Ti and the initial thermalization time, Tj for 
the background QGP expected to be formed at RHIC 
and LHC can be constrained to the total multiplicity as 
follows: 
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2tt 4 



45C(3) 4a eff 



dN 
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where Ra is the radius of the system, £(3) is the Rie- 
mann zeta function, a e ff = 7I ' 2 5e///90, and g e ff(= 
2x8 + 7x2x2x3x^/8) is the degeneracy of 
quarks and gluons in QGP and Np is the number of 
flavours. The value of the transition temperature, T c has 
been taken to be 175 MeV. We have used the boost in- 
variant model of relativistic hydrodynamics proposed by 
Bjorken [42| for the space time evolution of the expand- 
ing QGP back ground. The geometry of the collision and 
the space time evolution has been described in detail ear- 
lier, therefore, we avoid further discussions here and refer 
to our previous work [32j for details. The value of Ti and 
Tj for the QGP fireball are taken as T = 300 MeV and 
Ti = 0.5 fm/c. The corresponding quantities for LHC are 
Ti = 550 MeV and r t = 0.1 fm/c. The pressure (P)- 
energy density(e) relation for the QGP has been taken as 
P = e/3. 
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FIG. 2: Drag coefficients of bottom assuming running strong 
coupling, as(T) and temperature dependent Debye screen- 
ing mass, mo(T) due to its interaction with thermal gluons, 
quarks, and antiquarks. Inset: The variation of drag with p 
in the low domain. 



V. RESULTS 

With the formalism described and the inputs men- 
tioned above, we present the results now. The mo- 
mentum dependence of the drag coefficient (7) of charm 
quark propagating through QGP has been displayed in 
Fig. [1] for T = 300 MeV. The variation of 7 with p is 
non-negligible both for radiative as well as collisional pro- 
cesses. The drag coefficients, 7 is related to the energy 
loss as: 7 ~ (l/p)AE/Ax which in the high energy limit 
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(E ~ p) varies very slowly with E (Fig. [IJ , a behaviour 
very similar to the one observed in the fractional en- 
ergy loss as a function of the initial energy of light quark 
jets [43|]. The value of 7 due to collisional processes at 
p = 5 GeV is about 0.036 fm -1 which reduces to a value 
of 0.018 fm -1 at p = 10 GeV. Such a variation, which 
was neglected earlier, will have crucial consequences on 
the nuclear suppression factor, Raa for the charm and 
bottom quarks. In the inset of Fig. Q] we display the drag 
coefficient for collisional process for the low momentum 
domain. In the limited momentum range the drag re- 
mains almost independent of momenta. This constant 
values of drag has been used earlier in the FP equation 
and subsequently the solution was used in estimating the 
nuclear suppression factor. Since the value of 7 reduces 
with p one will overestimate the suppression by taking 
its value at low p. For the low momentum region we use 
7 = Jcoii- In Fig. [5] we depict the drag coefficients of a 
bottom quark which shows slower (compared to charm) 
variation with momentum. The drag coefficient for colli- 
sional loss has been shown in the inset of Fig. [2] for low 
momentum domain. 

The suppression of both charm and bottom quarks [32| 
(before fragmentation to hadrons) are plotted against pt 
in Figs. |3]and 0] respectively. We note that if one takes 
the drag to be momentum independent (or more precisely 
takes the value of 7 at low p and extends it upto very 
high p) then the drag due to collisional process causes 
about 50% suppression (dashed line). However, if we 
take into account the variation of 7 with p obtained from 
pQCD calculation then about 20% of suppression can be 
achieved, i.e. the contribution from the collisional loss 
becomes smaller with the momentum dependent drag. 
Therefore, the observed large suppression of the heavy 
quarks at RHIC is predominantly due to radiative loss. 
In fact, the inclusion of the radiative processes increases 
the suppression to about 75%. This can be understood 
from the fact that the drag due to the radiative loss is 
large. The suppression of the bottom quark is much less 
because of the smaller values of drag and initial harder 
momentum distribution (Fig. 2]). 

The hadronization of charm and bottom quarks to D 
and B mesons respectively are done by using Peterson 
fragmentation function [44J. The variation of the nu- 
clear suppression factor, Raa [12] with px of the D 
and B mesons has been displayed in Fig. [5] for RHIC 
initial condition (Tj = 300 MeV). The suppression for 
bottom is much less for the reasons mentioned earlier. 
The theoretical results show a slight upward trend for 
Pt above 10 GeV both for mesons containing charm and 
bottom quarks. Similar trend has recently been experi- 
mentally observed for light mesons at LHC energy [45| . 
This may originate from the fact that the drag (and hence 
the quenching) for charm and bottom quarks are less at 
higher momentum. 

The same formalism is extended to evaluate the nuclear 
suppression factor, Raa both for charm and bottom at 
LHC energy. Result has been compared with the recent 
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FIG. 3: Suppression of momentum of charm quarks in QGP 
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as a function px 



ALICE data(Ref. gg) in Fig.© The data is reproduced 
well by assuming formation of QGP at an initial temper- 
ature ~ 550 MeV after Pb+Pb collisions at v / sn"/v = 2-76 
TeV. 



VI. SUMMARY AND DISCUSSIONS 

The temperature and momentum dependence of drag 
and diffusion coefficients of heavy quarks interacting with 
the thermal partonic medium have been evaluated by us- 
ing the elastic and inelastic interactions. We have em- 
ployed both the dead-cone and the LPM effects in the 
calculation for the inelastic processes. The initial px dis- 
tributions for charm and bottom arc taken from MNR 
code [35| . The FP equation has been solved with the 
momentum dependent transport coefficients and subse- 
quently the nuclear suppression factors, Raa have been 
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FIG. 5: Raa as a function of pr for D and B mesons at 
RHIC. Experimental data taken from [3] and ||. 



calculated for D and B mesons for RHIC and LHC condi- 
tions. The momentum dependence of the drag coefficient 
is found to be crucial in reproducing the trend in the pr 
dependence of the experimental data. It has been seen 
that the radiative loss plays more dominant role than the 
collisional process. 
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